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DIFFERENTIAL FLATNESS AND COOPERATIVE TRACKING IN THE 

LORENZ SYSTEM* 

LUIS G. CRESPO* 


Abstract. In this paper the control of the Lorenz system for both stabilization and tracking problems 
is studied via feedback linearization and differential flatness. By using the Rayleigh number as the control 
-only variable physically tunable- a barrier in the controllability of the system is incidentally imposed. This is 
reflected in the appearance of a singularity in the state transformation. Composite controllers that overcome 
this difficulty are designed and evaluated. The transition through the manifold defined by such a singularity 
is achieved by inducing a chaotic response within a boundary layer that contains it. Outside this region, a 
conventional feedback nonlinear control is applied. In this fashion, the authority of the control is enlarged 
to the whole state space and the need for high control efforts is mitigated. In addition, the differentia 
parametrization of the problem is used to track nonlinear functions of one state variable (single tracking) 
as well as several state variables (cooperative tracking). Control tasks that lead to integrate and non- 
integrable differential equations for the nominal flat output in steady-state are considered. In particular, a 
novel numerical strategy to deal with the non-integrable case is proposed. Numerical results validate very 

well the control design. 

Key words, feedback linearization, differential flatness, Lorenz system, cooperative tracking, non-linear 
control 


Subject classification. Applied and Numerical Mathematics 
1. Introduction. The analysis and control of chaotic systems have attracted considerable attention 
in recent years. A chaotic system is strongly sensitive to small changes in the initial conditions. Such a 
behavior can be beneficial or detrimental depending upon the system and the objective under investigation. 
In this paper we will make use of the chaotic response of the Lorenz system to enhance the performance an 

controllability of feedback linearization based controllers. 

The Lorenz system is a simplified model of a thermally driven fluid convection system between parallel 
plates. Depending on the system parameters, such a system exhibits a rich spectrum of responses. The 
control of the Lorenz svstem has been studied by several researchers in recent years. Vincent and \u [ j 
proposed a bang-bang optimal control for stabilizing the unstable equilibrium points of the system^ Gao and 
colleagues [5] studied the nonlinear feedback control based on state space design. Chen and Liu [2|, lalwar 
and Namachchivava [9] and Alvarez-Gallegos [1] investigated the nonlinear regulation of the Lorenz system 
by using the feedback linearization techniques with different control structures and control objectives. 

From the physics point of view, it is natural to select the Rnyleigh number as the control variable. 
However, such a strategy makes the system uncontrollable at the plane * = 0. This control variable was 
used by the author in [3]. When feedback linearization techniques are applied, this feature appears as a 
singularity- in the state transformation. When the feedback linearization techniques are applied, this feature 
also appears as a singularity in the state transformation. The controllability of the system is limited to half 
of the state space. The singularity is also responsible for extremely high control efforts in its vicinity. In this 

*This~ work was supported bv the National Aeronautics and Space Administration under NASA Contract No. NAS1-97046 

while the author was in residence at 1CASE, NASA Langley Research Center lgcrespo0lc ase.edu. 
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paper, we propose composite controllers to overcome this limitation. Within a boundary layer that contains 
such a singularity, the chaotic response of the system and its response to step inputs are used to drive the 
system through. Furthermore, the differential flatness of the system is used to aim for control objectives 
that do and do not admit a closed form expression for the corresponding flat output of the system. 

This paper is organized as follows. A brief description and analysis of the Lorenz system are presented 
m Section 2. Section 3 discusses the feedback linearization and its application to the Lorenz system. In 
Section 4, composite controllers are designed and evaluated for both stabilization and tracking problems. 

2. The Lorenz System. The Lorenz model is obtained from studying a fluid layer heated from below 
and cooled from above such that a temperature difference is established across it. The convection motion 
is described by the Navier-Stokes equations. Taking Fourier expansion of these equations along two spa- 
tial directions and truncating the remaining expressions to retain only three modes leads to the following 
simplified model 


i = a{y-x), 

y = px — y — xz, (2.1) 

z = —dz + xy, 

where a, p and 3 are real parameters denoting the Prandtl number, the Rayleigh number and a geometric 
factor, respectively. The state variables x, y and 2 represent measures of fluid velocities and the spatial 
temperature distribution in the fluid layer under gravity. From the physical point of view, the Rayleigh 
number p can be easily manipulated by changing the heat transfer to the fluid from below. This parameter 
will be treated as the control variable. We denote u = p. 

For the brevity of discussion, we assume that the control converges to a constant value in steady state 
denoted as u ss = lim hoo u. The singular points of the system can be parametrized with u ss . The locus of 
these points is given by a family of curves 

^i = [0t0>0] T , X 23 = [±v/0(u„ - l),±\//?(u s ., - 1 ), Uss - 1] T . (2.2) 

Linearization about AT* leads to characteristic equation A 3 + .4A 3 + B\ + C = 0, where A = 1 -t- 3 + a 

B = 3{(J + 1) + a - u as and C = 3(<r-u as ). Linearization about X 2 3 leads to A = 1 + 3 + <t, B = 3(u aa + cr) 
and C = 23<j(u aa - 1). 

The stability analysis leads to the following observations. When u 88 < 1 the origin of the system is a 
stable equilibrium point. When 1 < u ss < p = a{o + 3 + 3)/(a -0-1), A',* is unstable and AT 2 * 3 are stable. 
When u ss > p there are no stable equilibrium points and the system reaches a chaotic regime. 

At Uss — 1 a pitchfork bifurcation takes place at AT* , while for u ss = p a subcritical Hopf bifurcation 
occurs at X^ 3 . For u ss > p, the system is driven by repulsions exclusively while the trajectories are confined 
to a region of finite volume forming a strange attractor. The response on the attractor is chaotic. For 

additional information, the reader can refer to [8]. In the numerical simulations, the parameters to be used 
are a - 10 and 3 = 8/3. 

3. Feedback State Linearization. 

3.1. Background. Consider the single input system 

X = f(x) + g(x)u, (3.1) 
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where x e R" is the state vector, u€ R is the control and f,g : R'-t R" are sufficiently smooth nonlinear 
functions of their arguments. The Lie derivative of g(x) with respect to the vector field f(x) is defined as 

arifg(x) = [f, adj -1 g](x) for k > 1 , 

a<F { g(x) = g(x), ( 3 - 2 ) 

where [X, Y] = VY ■ X-VX • Y is the Lie bracket of the vector fields X and Y. According to [6], there 
exists a real value function A(x), defined in a neighborhood U(x 0 ) of x 0 such that 

L g A(x) = L adtg X(x) = ... = L ad n-z g A(x) = 0, (3-3) 

L adr ‘g A ( X °) ^ °» 

where L g A(x) = VA(x) • g denotes the Lie derivative of the real-value function A(x) with respect to the 
vector field g, if (i) the matrix 

C = [g(xo), adfg(xo), adfg(xo), ..., ad 7 f 'g(x 0 )], (3-4) 

has rank n and (ii) the distribution D =span{g(x),ad f g(x),ad 3 g(x), g(x)} is involutive. Further- 

more. there exists a transformation z = $(x) in L r (x 0 ) such that 


<F(x) = [A(x), Lf A(x), ..., L? 1 A(x)] r , 

. x d n_1 A(x) T 
= [A(x), A(x), ] T - 


The system (3.1) is transformed into a chain of integrators 

Z\ = 2'2, Z'l — Z 3, • • • , Z n — ] -- Z n , Z n = fl(x) + b(x)lt = V , 


(3.6) 


n(x) = a(x) + b{x)u = L g L" *A(x) + If A(x)u, (3 7) 

where z* = Vf 1 A(x) (i = 1 n). The control for this linear system can be designed by full state feedback 

and pole placement techniques. For example, the control v can bo taken as. 

r($''(z)) = -Kz, (3.8) 

where the feedback gains K arc chosen to place the closed loop poles at the desired locations in the left hand 
side of the complex plane. For tracking problems, the feedback signal will be the tracking error z d -z where 
z d is a pre-specified reference signal in the z-domain. 

The control in the physical domain is then given by: 

u(x) = [t;($ _1 (z)) -«( x )] / b ( x )- 

It should be noted that u(x) becomes unbounded when 6(x) -t 0. 

3.2. Feedback Linearization of the Lorenz System. For the Lorenz system, f = [t r{y - x), -y - 

xz, -0z + xy] 7 and g = [0, x,0] T . After some manipulations we obtain 

adg(adfg) = [0,2<rx,0] / , 

arifg = [-a x,o{y - x) + x, ~z 2 ]' , (3.10) 

a (a - l)x - 2a 1 y 

udf g= -x 3 - 2axy + (a - l) 2 x + <r(l - a)y 

(3(7 - 1 - 3)x 2 - 2(7 xy 
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By evaluating the rank of the C matrix in Equation (3.4), we find that the rank is 3 except when the 
transformation is singular at x = 0 or a = 0/2. At the singularity at. x = 0,-the system (2.1) is completely 
insensitive to the control. This restriction imposes an unavoidable barrier in the controllability of the system, 
leaving just, half of the state space at the disposal of the control. Which half the system stays in depends on 
the initial condition. The closer the system gets to this singular plane, the higher the control effort will be, 
approaching positive or negative infinity as x approaches to zero. 

From Equation (3.3), the function A(x) , known as the flat output of the system [4, 7], satisfies 


OX 

X 7T = °> 
ay 


d\ 


ox 


<y)x + (T )-q-- x 


,ox 

dz 


= 0 . 


(3.11) 


Solving this set. of equations, we obtain A = x 2 /2 -a z + k, where k is the integration constant. Equation 
(3.5) leads to the state transformation 


x=[zi,z 2 ,z 3 ] T = $(x) (3.12) 

=[x 2 / 2 — az + k, u(0z — x 2 }, ajxy + 2a 2 x 2 — a0 2 z] T , 


where 7 = 0 — 2a. The inverse transformation is given by 


x=[a :,y,z}‘ = 4> l {z) 


4- Z 2 ) ^ 2a3z\ 4 (/? 4 2ct)z2 4 z% 


2a 3k 2a e 4 Z2 


n T 


os/2r00e + z^j 


aj 


where e — z\ - k. The transformed dynamic system takes the form of Equation (3.6) with n 
and 6(x) in Equation (3.7) given by 


(3.13) 


3 and a(x) 


a(x) = -ayx 2 z - 4a 3 x 2 - a { 7 4 ft 2 4 a (3 - 6a 2 )xy 4 a 2 'yy 2 4 a0 3 z , (3.14) 

b(x) = cryjr 2 . 


By using the solution of the PDE (3.11) and the system equation (2.1), we have 

A = x 2 /2 - az + k , (3.15) 

A = a(/3z — x 2 ), 

A = a'yxy 4 2a 2 x 2 - (3 2 az . 

The state variables and the control can be differentially parametrized using the flat output. 

X = 

y = ±(A + (0 + 2a)X + 2a 0\ - 2ak)/{ayj2~fx), (3.16) 

z = (0\ 4 2a 3\ - 2ak)/(a3j), 

U ~ 4a0yx + x 4 2 ca - ak) A - 23{Akc - 2 c 2 - 2 k 2 4 2 a^x 4 

,/?7X 4 2a d'ye - 7 a(3k)X - 4a3{3^x 4 4kc - 2 c 2 - 2k 2 )X - 8ake 2 4 4aj3~fk>c - 8ak 3 4 
I6ak 2 c 4 d 1 7(/J 4 2a)X 2 4 28 1 {3 4 <r)AA 4 .^A 2 4 2a/? 2 7 AA 4 2a/? 3 7 AA}, 

where c. = A 4 8X, x = c - k and v = A. Inequality range constrains on u can be imposed by designing 
composite controllers. This is explained in the next section. 
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4. Composite Control. From the above discussions, we know that a feedback linearization based 
control is not able to drive the system across the singularity imposed by the transformation. This fact 
restricts the controllability of the system to a portion of the state space, half in this case, and leads to 
extremely high control effort in the vicinity of the singularity. In i his section we propose a composite 
controller to overcome such a difficulty. 

Denote the hyper-plane that makes the state transformation singular as g(x) = 0. Define the composite 
controller given by: 


u\ if |x - < 7 ( x ) j > 6 
u 2 otherwise 


(4.1) 


where 6 defines the thickness of a boundary layer about the singular plane, u x is the control expression 
resulting from using a conventional control method and u 2 is the amplitude of a step input set according to 
the particular control objectives. 

The condition for applying u\ can be replaced by |ui| < L\ provided that control saturation occurs in 
the vicinity of the boundary layer. In this paper, g(x) = 0 and u 2 is a step input that induces a chaotic 
response within the boundary layer i.e. u 2 > p . 

The non-empty intersection of the attractor arid the two hyper-planes defined by |x-<?(x)| = S guarantees 
that the crossing of the boundary layer occurs. This can be proved as follows. Let’s call hi(x) for i = 1,2 
these two planes. Assume that (i) the chaotic response is moving within a strange attractor whose state 
space location is given by A\ and that (ii) ,4i n /i,-(x) 0. The crossing of the hyper-plane /i*(x) will not 

occur iff there exist a subset A 2 C A\ such that A 2 Pi /?,(x) — 0 for all times i.e. A 2 is a strange attractor 
by itself. This implication violates the irreducibility property of .4] i hen A 2 can not exist. The need for 
A\ fl hi(x) ^ 0 imposes bounds to S from above. 

Notice that in this scheme the control does not have the authority to manipulate the transient part of 
the transition from one side of the boundary layer to the other one. For some states several crossings of g{x) 
might occur before the system leaves the boundary layer. Such behavior is clearly undesirable. Refinements 
and improvements of the control within the boundary layer can be achieved by taking into consideration 
the control objective and the state of the system. This practice however was not implemented in this study. 
Once the system leaves the boundary layer ii\ is applied and the stabilization/tracking is achieved. In this 
paper u\ is given by u in Equation (3.16). Due to the structure of the controller, global uniform asymptotic 
stability about the reference signal r(x(f)) is achieved at the locations where the intersection of r(x(t)) and 
the boundary layer is an empty set. 

For stabilization, the linear system given by Equations (3.6), (3.7) and (3.14) can be controlled by pole 
placement techniques. Taking the feedback control v as: 


v = a\Z\ + a 2 z 2 4- o 3 z 3 , 


(4.2) 


where the feedback gains a j are chosen to place the closed loop poles in trie left hand side of the complex plane. 
On the original state variables, the control can be obtained after substii uting Equations (3.12) and (3.8) into 
Equation (3.7). Since this procedure stabilizes z, the steady-state values in the ^-domain can be controlled by 
manipulating k according to Equation (3.13). The reader must notice that stabilization about the origin using 
feedback linearization requires infinite control effort, i.e. lim^o u(x) — lim T _>o(^(x) — a(x))/b(x) = ± 00 . 
However, from the stability analysis we know that any control satisfying 0 < u S8 < 1 will drive the system 
to the origin from a given initial condition. 
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4.1. Single State Tracking Control. In this section we present examples for stabilization and track- 
ing control problems that involve a single state. The tracking signal x d (t,) = a + bsin(t) for the state variable 
x(t) is considered here. The flat output corresponding to the tracking signal x d {t) is given by 

A d {t) = h/2/3)(b 2 /2 + a 2 ) + ('yab/(P 2 + *)) {/?sin(<) - cos (f)} (4.3) 

- (b 2 y/(2(3 2 + 8)){(/3/2) cos(2t) - sin(2i)} 

The difference between the system flat output, X d - A = z d - z\ is the tracking error. We take the full state 
feedback control v for the tracking problem as follows 

v = Xd + A i (Arf — A) + K 2 (X d — A) + A 3 (A^ — A) (4.4) 

The control gains Kj arc selected such that the tracking error vanishes exponentially. It should be noted 
that the tracking control is designed in the transformed space z, and therefore is indirect for x. 

The state variables and the nominal control obtained from Equation (3.16) are shown in Figure 5.1 for 
a = 5 and 5 = 6. The control grows unbounded near the singularity. Time evolutions for different x d (t) 
are shown in Figures 5.2 and 5.3 using S = 0.1 and 6 = 0.2 respectively. In all the cases the control is 
activated after 30 seconds of chaotic regime. In the first case, the system does not reach the boundary layer 
and perfect tracking is achieved after a short transient. If the control is activated when the singular plane 
is between the state of the system and x d (t), the system would reach and cross the boundary layer before 
settling down. This control was designed such that u 2 is applied when u\ < 6. In this fashion the control 
range constraint u > 0 is imposed. In the second case, the desired trajectory crosses repeatedly the singular 
plane forcing the system to reach the boundary layer several times. The effect of not applying the nominal 
control is slightly noticeable. The reader must notice that both the x and y states behave similarly once the 
system reaches the strange attractor. 

Figure 5.4 shows the results of the composite controller with a wider boundary layer S = 0.5. Recall 
that the control does not have authority within the boundary layer and relies on the chaotic behavior of 
the system to cross it. As expected, the increase of 6 has a detrimental effect on the control performance. 
For even higher values of 5, the system might be trapped in the boundary layer while the desired trajectory 
completes half a cycle. 

4.2. Cooperative Tracking Control. Now we use the differential parametrization of the states and 
control to aim for tracking objectives that involve combinations of the states. The problem statement is as 
follows. Find u\ in Equation (4.1) such that the system is driven to the manifold defined by h(xM,t) = 0 
from any initial condition. In this problem, the system would track signals that imply cooperative relations 
among the states, being the tracking of a trajectory of a particular state a particular case. Again, the 
composite structure of the control enable to achieve global stability about the desired tracking function. In 
the examples to come we will take u 2 = 30 > p within the boundary layer. 

Depending upon the tracking objective, the equation for the desired flat output can be integrable or 
non-integrable. Both cases are considered next. 

4.2.1. Problems with closed form solution for X d . Once the control objective /?(x,u,£) = 0 is set, 
the Equations (3.16) lead to an ODE for the desired flat output A d . Such an equation is in general non-linear. 
In this section we study problems where a closed form expression for the steady state solution can be found. 

As the first example, we take h(x,u,t) - x 1 ~ z - a - bsin(ujt). The corresponding differential equation 
for the flat output is given by 

A d T (2<r£/0)A rf = 4- (fory/^) sin(u;t), 
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(4.5) 


where 0 = 2<r - 1 and £ = 9 - 1- Solving we find 

A d = 07/(20 + ?;exp{-2t<7£/0} + {26<t 2 7£ sinM) - uiba-fC cosM)}/(4ct 2 £ 2 + ^ 2 0 2 ) 


For our purposes, providing that the transient response vanishes, only the steady state component is needed. 
This expression along with Equations (3.16) and (4.4) fully determine the control. After some manipulations 
we find that the system reaches stationarity at 


x ss 


T x/tfa/O T \//la/0 a/^j . 
2 + (a — 0/£t 


(4-6) 


Numerical results for a — 
the boundary layer twice. 


15 and 6 = 5 are shown in Figure 5.5. For this particular case, the system crosses 
Figure 5.6 shows the numerical results for a = 10 and 6 = 0. 


As a second case, we use h(x) = * 2 + .V 2 + * 2 - This problem can be interpreted as the stabilization 
of the system about the surface of a sphere. The corresponding differential equation for A d , not shown here, 
is non-linear. However, because the homogeneous solution vanishes with time we can solve for particular 
solution and use it in Equation (4.4) to calculate the control. The corresponding flat output and steady 

state values are: 


Xd = 7(±# - /3)/2, 


X g s — 


±y/ 0 (O- 0),±\/0 (0 - 0 ), 6 - B ] > 


».,■ = ± (1 + 0 - 9 ) , 


where 6 = s/9 2 + /? 2 . Stationary values of the states can also be found by solving for the values of k in 
Equation (3.13) that satisfy h(x) = 0 when z -+ 0 (due to Equation (3.8). Numerical results for /? = 2 are 
shown in Figure (5.7). In the case shown, the boundary layer is not. reached. 

4.2.2. Problems without closed form solution for X d . In the previous section, the tracking objec- 
tives led to stable ODEs for A d whose particular solution could be found in closed form. In this section, we 
consider problems in which this is not the case. 

By integrating numerically and simultaneously both the state Equations (2.1) and the ordinary differen- 
tial equation for the desired flat output A d , tracking can be achieved. It is important to notice that realizable 
objectives imply stable solutions for A d . The tracking error dynamics, set by Equation (4.4), makes the flat 
output in Equation (3.15) converge to the steady state value of A d . Tracking is achieved once the transient 
for both the real and the desired flat outputs vanish. 

Due to the non-linear character of the equation for A d , the steady state response might exhibit dependence 
to the initial conditions. For this reason, the control design must start by searching for the initial conditions 
in A d and its derivatives that lead to steady state trajectories that saiisfy the desired response specifications. 
The reader must notice that a fixed number of derivatives of A d , three in this case, are needed to build ». 
This can be obtained by performing additional differentiations on the ODE for A d . 

As an example, we use the energy-like expression h(x,t) = + V 2 + mz - E(t). The corresponding 
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differential equation for A,/ is given by: 


0 = -2E3y 2 a 2 (X d + A d 0 - k) + 2mya3X' 2 d - 4mkya 2 X d + 4mya 2 0X d X d + 

4mya 2 0 2 X d + 4mk 2 ya 2 + 2mya3 2 X d X d - 2myafikX d - 8mya 2 3kX d + 

24 fi 2 k 2 a 2 X d - 24a 2 0kX d + 24A : 2 <j 2 3X d - 48a 2 k3 2 X d X d + 8a 2 3X d - (4.7) 

8k 2 <j 2 3 + 43 3 ya 2 X 2 d + y3X 2 d + 2y3 2 X d X d + 4ya3X d X d + 4ya0 3 X d X d - 
8kya 2 3X d + 4ya0 2 X d + 4k 2 y a 2 3 - 4kya3X d - 4kya3 2 X d + 8ya 2 0 2 X d X d + 

4ya 2 0X 2 d - 83 2 kya 2 X d + 4ya3 2 X d X d + y3 2 X d + 83* o 2 X 2 d + 248 3 a 2 X 2 d X d - 
243 A a 2 kX 2 d + 24a 2 3 2 X d X d 

The set of equations to be integrated numerically is given by Equations (2.1) and (4.7), where u is given 
by Equations (4.1) and (4.4). The reader must notice that does not depend on the states and that both 

equations are coupled via u. Numerical simulations with E(t) = a + 6sin(wt) cos(«,<) are shown in Figure 
(5.8). 


5. Conclusions. This paper studies the stabilization and tracking control of the Lorenz system using 
feedback linearization and differential flatness. When the Rayleigh number is used as the control variable, 
the system is uncontrollable in a manifold of the state space. In the vicinity of such a singularity, the control 
demands grow unbounded. Composite controllers that use feedback linearization and the svstem response 
to step inputs are proposed to overcome this difficulty. By inducing the chaotic response within a boundary 
layer which contains the singular plane, the transition to desired states is achieved. Such controls can be 
used not only to enlarge the controllability region of the system to the whole state space but also to mitigate 
high control demands. Control objectives and initial conditions that imply single and multiple crossings of 
the boundary layer are studied in the examples. In addition, tracking control problems that involve single 
and cooperative relations among the states are studied using the differential flatness of the system. Problems 
with control objectives that lead to integrable and non-integrable differential equations for the desired flat 
output are considered. A numerical approach in which the state equations and the differential equation for 
the nominal flat output are simultaneously integrated is proposed and validated. Numerical simulations led 
to excellent performances. 
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Time 

Fig. 5.3. Time evolutions of the state variables and the control with 6 = 0.3, a = 0 and b = 10. The control is activated 
after 30s. Sharp peaks in the control curve indicate that the system reaches the singularity . 



Time 


Fig. 5,1. Time evolutions of the state variables and control with S = 0.5, a = 0 and b = 10. The control is activated after 
30s. Sharp peaks in the control curve indicate that the system reaches the singularity. 
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100 



Fig. 5.5. Time evolutions of the state variables and the control for tracking the hyper-plane h(x , t) = x 2 -z- 15- 
The following conventions are used xfdotted line ), y(dot-dashed line ), z(dashed line), ufsolid line) and hfthick line) 



Fig. 5.6. Time evolutions of the state variables and the control for tracking the hyper-plane /i(x) = x 2 — z — 
following conventions are used xfdotted line), y( dot- dashed line), z(dashed line), u(solid line) and h(thick line). 
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5sin(u;£). 


10. The 
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Fic 5 7 Time evolutions of the state variables and the control for tracking the hyper-plane ft(x) - x 2 + y 2 + z 2 4. The 
following conventions are used x (dotted line), y(dot-dashed line), z(dashed line). u(solid line), and hfthick line) 



Fig. 5.8. On the top , time evolutions of the states x (dotted line), 
middle , time evolutions for the flat output and its derivatives are shown, 
desired (solid line) and the real (dotted line ) flat outputs are shown . 


y (dashed line) and z (solid line) are shown . In the 
In thf bottom, time evolutions of E(t) for both , the 
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